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Summary 

This study will demonstrate the importance of application of CNC machines in generation 
of gear tooth surfaces with new topology. This topology decreases gear vibration and will 
extend the gear capacity and service life. A preliminary investigation by a tooth contact 
analysis(TCA) program has shown that gear tooth surfaces in line contact (for instance, 
involute helical gears with parallel axes, worm-gear drives with cylindrical worms etc.) are 
very sensitive to angular errors of misalignment that cause edge contact and an unfavorable 
shape of transmission errors and vibration. 

The new topology of gear tooth surfaces is based on the localization of bearing contact, 
and the synthesis of a predesigned parabolic function of transmission errors that is able to 
absorb a piecewise linear function of transmission errors caused by gear misalignment. 

The report will describe the following topics: (1) Description of kinematics of CNC 
machines with 6 degrees-of-freedom that can be applied for generation of gear tooth surfaces 
with new topology. (2) A new method for grinding of gear tooth surfaces by a cone surface 
or surface of revolution based on application of CNC machines. This method provides an 



optimal approximation of the ground surface to the given one. This method is especially 
beneficial when undeveloped ruled surfaces are to be ground. (3) Execution of motions of 
the CNC machine. The solution to this problem can be applied as well for the transfer of 
machine-tool settings from a conventional generator to the CNC machine. 

The developed theory required the derivation of a modified equation of meshing based 
on application of the concept of space curves, space curves represented on surfaces, geodesic 
curvature, surface torsion etc. Condensed information on these topics of differential geometry 
is provided as well. 


Introduction 

The design and manufacture of gears with new topology of gear tooth surfaces are prob- 
lems of great importance for helicopter transmissions. The existing technology of gears is 
restricted with the necessity to use cutting and grinding machines whose kinematics is based 
on linear relations between the motions of the tool and the workpiece. 

The need of low-noise gears with increased load capacity and service life can be satisfied 
with a new topology of gear tooth surfaces that is able to provide: (i) a reduced sensitivity to 
misalignment and avoidance of edge contact, (ii) a parabolic type of function of transmission 
errors to reduce the level of possible vibration, (iii) a localized bearing contact with controlled 
dimensions of the instantaneous contact ellipse. 

The application of CNC (Computer Numerically Controlled) machines overcomes the 
obstacles presented for generation of gears with a new surface topology by using the existing 
equipment. The CNC machines are able to provide computer controlled nonlinear relations 
between the motions of the tool and the gear being generated. Although such machines are 
used at present mainly for the installment of machine-tool settings with higher precision, 
their prosperous future is in their application for generation of gears with new topology. The 
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CNC machines are a unique opportunity for researchers to modify the geometry of traditional 
gear drives and benefit industry with gear drives with substantiaUy improved parameters. 

The modification of geometry of gear tooth surfaces requires from researchers a new 
approach for the development of principles of conjugation of gear tooth surfaces. Application 
of conjugate gear tooth surfaces being in instantaneous line contact is in the authors’ opinion 
an anachronism. Such gear tooth surfaces are very sensitive to misalignment that causes the 
shift of the bearing contact to the edge and transmission errors of such a type that cause 
a jerk at the transfer from one cycle of meshing to the next. The considerations above are 
true for spur gears, involute helical gears with parallel axes, and worm-gear drives. 

It is necessary as well to change the attitude to some finishing processes such as honing 
and shaving applied for helical and spur gears. It is not reasonable to require that such 
finishing processes would provide the exact screw involute surfaces knowing ahead that only 

modified tooth surfaces are to be applied. 

The statements mentioned above are illustrated with the following drawings. 

Figures 1 to 4 show the influence of angular errors of misalignment A 7 and d\ pi of involute 
helical gears with parallel axes that cause edge contact (figures 1 and 3) and piecewise almost 
linear functions of transmission errors (figures 2 and 4). The drawings above are based on 
the investigation performed by Reference [1]. The design parameters of the helical gears are 
shown below: 
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Figures 5, 6 and 7 illustrate the impact of misalignment of worm-gear drives. The shift 
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of the bearing contact is shown in figures 5 and 6, and the undesirable shift of transmission 
errors is shown in figure 7 that will inevitable cause vibration. The drawings are based on 
the research accomplished by Reference [2]. 

The design parameters of the worm-gear drive are shown below: 


Nr 

= 2 

n 2 

= 30 

Axial module, m 

= 8(mm) 
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II 
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Shortest distance, E, 

= 176(mm) 


Figure 8 illustrates why a predesigned parabolic type of transmission errors is beneficial 
for the gears with the new topology (Ref. [3]). This figure illustrates the interaction of 
a parabolic function of transmission errors (provided at the stage of synthesis of the gear 
tooth surfaces) with a linear function of transmission errors (caused by misalignment). The 
combination of these functions is again a parabolic function, with the same slope as the 
predesigned one that is translated with respect to the initial parabolic function. This means 
that the predesigned parabolic function absorbs the linear function and keeps the shape of 
a parabolic function. 

There are three cases of generation of the workpiece surface £ p by the given tool surface 
£t by CNC machines: 

(1) Surfaces £ t and £ p are in continuous tangency, however they contact each other at 
every instant at a point not a line. 

(2) Surfaces £ t and E p are in continuous tangency and they contact each other at every 
instant at a line. Surface £ p is generated in this case as the envelope to the family of 
surfaces E t . The family of surfaces is generated in relative motion of E t to E p . 

(3) An approximate method for generation of a surface T, g (ground or cut) with an optimal 
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approximation to the ideal surface E p . 

An example of case 1 is the generation, for instance, of a die designed for forging of a gear. 
Generation of conventional spiral bevel gears and hypoid gears by the Phoenix machine 
is the example of case 2 generation. Case 3 is the basic idea for a new method for surface 
generation discussed in section 4. Only cases 2 and 3 of surface generation are discussed in 
this report. 

The contents of the report covers the following topics: 

(i) Description of “Phoenix” and “Star” CNC machines, that are suitable for generation 
of gear tooth surfaces with new topology. 

(ii) Execution of motions of CNC machines. 

(iii) Generation of a surface with optimal approximation to the ideal surface. 

(iv) Concept of curvatures that are required for computations for the proposed approach 
for generation. 


1. “Phoenix” and “Star” CNC Machines 
“Phoenix” CNC Machine 

The “Phoenix” CNC machine (figure 9) is designed by the Gleason Works for generation 
of spiral bevel and hypoid gears. The machine is provided with a total of six degrees-of- 
freedom. Three rotational motions, and three translational motions are used. The transla- 
tional motions are performed in three mutually perpendicular directions. Two of rotational 
motions are provided as rotation of the workpiece and the rotation that enables to change 
the angle between the axes of the workpiece and the tool. The sixth rotational motion is 
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provided as rotation of the tool about its axis, and generally is not related with the pro- 
cess for generation. The motions with other five degrees-of-freedom are provided as related 
motions in the process for surface generation. 

Coordinate Systems Applied for “Phoenix” 

Coordinate systems S t (x t , y t , z t ) and S p (x p , y p , z p ) are rigidly connected to the tool 
and the workpiece, respectively (figure 10). For further discussions we will distinguish four 
reference frames designated in figure 9 as I, II, III and IV. The reference frame IV is the 
fixed one to the housing of the machine. Reference frames /, II and III perform translations 
in three mutually perpendicular directions, respectively. We designate coordinate systems Sh 
and S m that represent reference frames I and III , respectively(figures 9 and 10). Coordinate 
axes of Sh and S m are parallel to each other and the location of Sh with respect to S m is 
represented by (x^ h) , y^ h \ and Coordinate system S t performs rotational motion 

with respect to Sh about the z^-axis. To describe the coordinate transformation from S m to 
S p , we use coordinate systems S e and 5^ (figure 10). Coordinate system S e performs rotation 
with respect to S m about the y m -axis. Coordinate axes of system Sd are parallel to the 
respective axes of the location of origin Od with respect to O e is determined with the 
parameter x^°^ = const. Coordinate system S p performs rotational motion with respect to 
Sd about the x^-axis. 

“Star” CNC Machine 

A version of the “Star” CNC machine that is provided with 6 degrees-of-freedom is shown 
in figure 11. Coordinate systems S t (x t , y t , z t ), S p (x p , y p , z p ) and Sj (x /5 y f , z f ) are rigidly 
connected to the tool, workpiece and frame, respectively. Coordinate system Sd is parallel 
to system 5/ and the location of Sd with respect to 5/ is represented in . Sf by {x^ d \ 0, 
0). Coordinate system S e performs rotational motion with respect to Sd about the y^- axis. 
Coordinate system Sh is parallel to S e and the location of Sh with respect to S e is represented 
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in S e by (0, y[° h \ z£ 0h) )- Coordinate system S t performs rotational motion with respect to 
S h about the x*-axis. Coordinate system S p performs rotational motion with respect to 
the fixed coordinate system Sj about the x/-axis. Altogether there are three translational 
motions along axes xy , y e and z t and three rotational motions about axes xy, y^, and x^. 

2. Basic Principle of Execution of Motions on CNC Machine 

Consider that the location and orientation of the tool with respect to the workpiece are 
given in coordinate systems that are represented for a conventional generator or for an ab- 
stract (mathematical) model of the process for generation. We will consider for the following 
derivations the example of application of the “Phoenix” machine. A similar approach can be 
applied for other types of CNC machines, for instance, for the “Star” machine. Our goal is 
to develop the algorithm for the execution of motions of the CNC machine using the initial 
information mentioned above. Reference [4] has used for this purpose the existence of a 
common trihedron for the two couples of coordinate systems ( St , Sj, and (St » S p ) 
that are applied for the CNC machine and for the generating process, respectively. The 
approach used in. this report is as follows: 

(i) Consider that 4x4 matrices M^ and 3x3 matrices (k = C, G ) have been 
derived. The superscripts “C” and “G” indicate the CNC machine and the abstract 
generating process, respectively. 

(ii) The matrix equality 




( 1 ) 


will provide the same orientation of S ( t k) with respect to S^, k) ( k = C, G ) in both 
reference frames. 
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(iii) The matrix equality 


Mj?[ 0 0 0 1 ] r = Mj?[0 0 0 1 ] T 


( 2 ) 


will provide the same position vector ( 0 p 0 t ) p for both reference frames. 

The application of equations (1) and (2) for the execution of motions of the “Phoenix” 
machine is considered for the two following cases: (i) a hypoid pinion is generated by appli- 
cation of a conventional generator, and (ii) a surface with optimal approximation to the 
ideal surface E p is generated. 

Derivation of Matrix L p ^ and Position Vector ( O t O p 

Using a routine procedure for coordinate transformations, we obtain 


L — L p< i(0)L(ieL em (^)L TO ^L/ l i(/i) 


cos /j. cos 4> — sin \i cos <p sin <f> 

cos /x sin <{> sin xjj sin fi sin <j> sin t b cos 6 sin xjj 

+ sin fi cos i/> -f cos /x cos t /> 

cos /x sin <j> cos 0 sin /x sin <j> cos */> cos q cos 

— sin /x sin ip — cos ji sin t/> 


( 3 ) 


We note that Lj e and L m h are unit matrices. 

The derivation of the position vector (O t O p )^ in S p is based on the following consider- 
ations: 
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(0 


{O m O t )f' 1 + (0 t 0 p )< c) = (O m O p )J, c) 


Thus: 


(0^)< c > = (0ZU~ r )l c) 1 - (0^J7)< C > = (O.O,)' 17 ' - (O m O*)< c > 




(4) 


Here: x { ° d) = const., x [ ° h \ y ( ° h] and z£p h) are considered as algebraic values. 

(ii) Vector (O t O P ) (C) can be represented in coordinate system S* c) with the following matrix 

equation 


[O t O v )f ) = x[ 0i) i p - 1 0 Of 


- y ^L pm [0 1 Of-^VtO 0 if 


(5) 


where L pm = L pi L ie L ein (L* is a unitary matrix). 
Equation (5) yields 


( o,o r )l C) = 


x(°<d — x^ h ^ cos <p — z^f h ) sin <p 
sin <p sin ip - y£ h) cos ip - z(? h) cos <p sin ip 
x { ° h) sin 6 cos ip + sin ip - z<° h) cos <p cos ip 


( 6 ) 
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3. Example: Generation of Hypoid Pinion by “Phoenix” 


Generation of Pinion Tooth Surface by Conventional Generator 

The pinion tooth surface is generated a s the envelope to the family of tool surfaces that 
are cone surfaces as shown in figure 12. 

Henceforth, we will consider the following coordinate systems: (i) the fixed ones, S 0 and 
S q that are rigidly connected to the cutting machine (figures 13 and 14); (ii) the movable 
coordinate systems S c and S p that are rigidly connected to the cradle of cutting machine 
and the pinion, respectively; (iii) coordinate system S t that is rigidly connected to the head 
cutter. In the process of generation the cradle with S c performs rotational motion about the 
z 0 -axis with angular velocity and the pinion with S p performs rotational motion about 

the x,-axis with angular velocity u>^ (figure 14). 

The tool (head- cutter) is mounted on the cradle and performs rotational motion with the 
cradle. Coordinate system S t is rigidly connected to the cradle. To describe the installment 
of the tool with respect to the cradle we use coordinate system Sb (figures 12 and 13). The 
required orientation of the head-cutter with respect to the cradle is accomplished as follows: 
(i) coordinate systems Sb and S t are rigidly connected and then they are turned as one rigid 
body about the 2 c -axis through the swivel angle j = 2x - <5 (figure 13); (ii) then the head- 
cutter with coordinate system S t is tilted about the yj,-axis under the angle i (figure 12(b)). 
The head-cutter is rotated about its axis z t but the angular velocity in this motion is not 
related with the generation process and depends only on the desired velocity of cutting. 

The pinion setting parameters are E m - the machine offset, 7 m - the machine-root angle, 
A B- the sliding base, and A A- the machine center to back are shown in figure 14. The 
head-cutter settings parameters are Sr- radial setting, 6 C - initial value of cradle angle, j- 
the swivel angle(figure 13), and i- the tilt angle(figure 12(b)). 

Pinion Tool Surface Equations 
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The head-cutter surface is a cone and is represented in S t (figure 12(a)) as 

(r c + s sin a) cos 0 
(r c -f ssin a) sin0 

r,M) = 

—s cos a 
1 

Here: (s, 9) are the Gaussian coordinates, a is the blade angle and r c is the cutter point 
radius. Vector function (7) with a positive and a negative represents surfaces of two head- 
cutters that are used to cut the pinion concave side and convex side, respectively. 

The unit normal to the head-cutter surface is represented in S t by the equations 

n t = [ — cos a cos 0 — cos a sin 0 — sina] T (8) 

The family of tool surfaces is represented in S p by the matrix equation 

r p (s, 9, <f> p ) = M p , M gn M n0 M oc M c & M** r t (s,0) (9) 

Here: S n is an auxiliary fixed coordinate system whose axes parallel to S 0 axes. 

cos i 0 sin i 0 

0 10 0 

M bt = 

— sin i 0 cos i 0 

0 0 0 1 
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— sinj —cos j 0 Sr 

cos j — sin j 0 0 

M c6 = 

0 0 10 

0 0 0 1 

cos q sin q 0 0 

— sin q cos q 0 0 

M oc — 

0 0 10 

0 0 0 1 

10 0 0 

0 10 E m 

0 0 1 -A B 

0 0 0 1 

cos 7 „, 0 sin 7 m -A A 

0 10 0 

— sin 7 m 0 cos 7 m 0 
0 0 0 1 

10 0 0 
0 cos 0 V — sin 4> p 0 

M pg = 

0 sin (f> p cos 4> p 0 

0 0 0 1 

6 = 2w — j\ q = 9 e + rn C p4> v where 9 C is the initial cradle angle and rricp = /a/ p ). 

Equation of Meshing 

This equation is described in Reference [5] as: 




n (p) . v ('p) = N (p ) • v (cp) = /(s,M P ) = 0 


( 10 ) 


12 



where stud are the unit normal and the normal to the tool surface, and ^ is the 
relative velocity between the tool surface and the workpiece. 

Equation (10) is invariant with respect to the coordinate system where the vectors of the 
scalar product are represented. These vectors in our derivations have been represented in S 0 
as follows, 


n 0 = L oc L c b Lbt n t 

= [( u[ c) - x r 0 ] + (0^4 x w< p >) 


Here: 


r 0 = M oc M cb M fct r t 


O^A = [ 0 - E m A B] t 

= — [ cos 7 0 sin7] T ; (I w i p) | = 1) 

u>i e) = -[ 0 0 rricp ] T 

Pinion Tooth Surface 

Equations (9) and (10) represent the pinion tooth surface in three-parametric form with 
parameters s, 6 and 4> p . However, since equation (10) is linear with respect to s we can 
eliminate s and represent the pinion tooth surface in two-parametTic form as 

r„(0, i,,dk) (11) 
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Here: d k (k = 1, . . . , 8) designate the installment parameters: E m , 7 m , A B, A A, Sr, O c , 
j and i. 

The normal to the pinion tooth surface is represented as 

n p { 0 , 4 > p ,d k ) (12) 

where d k (k = 1,2, 3 , 4 ) designate the installment parameters 7 m , 6 C , j and i. 

Derivation of and (0 t 0 p ) j, G) 

Our next goal is to derive the algorithm for execution of motions on “Phoenix ”, knowing 
the basic machine-tool settings on the conventional generator. 

The coordinate systems applied for the CNC machine are represented in figure 10 . The 
performed coordinate transformation yields: 

(L pf ) (G) = M<z)] (k = 1 , 2 , 3 ;/ = 1 , 2 , 3 ) ( 13 ) 

Here: 

a n = cos * cos 7 m sin(<? — j) — sin i sin7 m , 

a 12 = ~ cos (q - j) cos 7 m , 

a 13 = sin i cos 7 m sin(^ — j) + cos isin7 m , 

U21 = cos i sin 7 m sin <t> p sin (q — j) + cos i cos (q — j) cos <f> p 
+ sin i cos 7 m sin 4 > p , 

a 2 2 = - cos (q - j) sin7 m sin 4 > p + sin (q - j ) cos (f> p , 

a 2 3 = sinz sin7 m sin <ji>p sin(<? — j) + sini cos(<7 — j) cos <^>p / ( 14 ) 

— cos i cos 7 m sin <(> p , 

CZ31 = — cos i sin 7 m cos <f> p sin(g — j ) + cos i cos(q — j) sin <j> p 

— sin i cos 7 m cos <j> p , 

a 32 = sin 7 m cos<f> p cos(q - j) + sin (q - j) sm<f> p , 
a 33 = — sin i sin 7 m sin(q — j) cos d p + sini cos(^ — j) sin <f> p 
+ cos 1 cos 7 m cos <f> p 
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The variable parameters q and <j> P are related and therefore the coefficients a kl ( k = 
1,2,3;/ = 1,2,3) are functions of q. 

The position vector (O t O p )jf is represented as follows: 


mf = -(M pl ) (G) [ 0 0 0 if 

Sr cos q cos 7 m — A B sin 7 m — A A 
— 5 fl(sin q cos <j> p — cos q sin 7 m sin <j> p ) 
+E m cos <f> P + Aj B cos 7 m sin 4> p 
— Sr( sin q sin <f> p T cos q sin 7 m cos <f> p ) 
+E m sin <j) p - A B cos 7 m cos <f> p 
1 

a 14(?) 

a 24(9) 

<134(9) 

1 


(15) 


Execution of Motions of CNC 

Matrix equality (equation (1)) provides nine dependent equations for determination of 
functions <£(<?), ?/>( q ), and fi(q). We can determine these functions using the following proce- 
dure: 

Step 1: Determination of 4>. 


sin <f> = a l3 (q) 

This equation provides two solutions <j > ; the smaller value of <j> can be chosen. 


( 16 ) 


15 



Step 2 : Determination of ifr. 


cos <j> sin — a 23 (g) , cos cos = <133(9) (17) 

These equations provide a unique solution for yj, considering <p as given. 

Step 3: Determination of y. 

cos y cos <j> = o n (9) , — sin ^ cos <^ = 012(9) (18) 

These equations provide a unique solution for y, considering <f> as given. 

For the generation a face-milled hypoid pinion, a tool with a cone surface is applied. The 
tool surface is a surface of revolution and the rotation of the tool about its axis is not related 
with <f>. Functions (17) must be applied and executed only for the generation of face-hobbed 
hypoid pinion, that is cut by a blade. 

Vector equality 


(0,0 r )W = (OtO r )P (19) 

permits the determination of functions x^ h \q), yl^ h \q), and zj^^(q). Equations (6), (15) 
and (19) considered simultaneously, represent a system of three linear equations in the un- 
knowns: The solution to these equations enables to determine the trans- 

lational motions on the CNC machine. 

4. Generation of a Surface with Optimal Approximation To the Ideal Surface 
Introduction 

This section is based on the research accomplished by Reference [2] that was directed at 
generation of a surface (E a ) that must be in optimal approximation to the theoretical (ideal) 
surface E p . 
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The method for generation of E a is based on following ideas: 

(1) A mean line L m on the ideal surface £ p is chosen as shown in figure 15. 

(2) The tool surface £ t is a properly designed surface of revolution (in particular cases £ t 
is a circular cone as shown in figure 15) that moves along L m . Surfaces St and £ p are 
in continuous tangency along L m \ M is the current point of tangency (figure 15). The 
orientation of E t with respect to £ p (determined with angle 0) is continuously varying. 
Angle 0 at current point M of tangency is formed by the tangents t / and t b to L m and 
the tool generatrix, respectively (figure 15). Tangents tj and t(, form plane II that is 
tangent to St and S p at point M. 

(3) The tool surface E t in its motion with respect to E p swept out a region of space as a 
family of surfaces E t . The envelope to the family of E t is surface E a , the ground or 
cut surface, that is in tangency with the theoretical surface E p at any point M of L m 
and must be in optimal approximation to E p in any direction that differs from L m - 

(4) The optimal approximation of E a to E p is obtained by variation of angle 0 (figure 15). 

(5) The continuous tangency of tool surface E, with E p and properly varied orientation 
of E t can be obtained by the execution of required motions of the tool by a computer 
controlled multi-degree-of-freedom machine. One of these degrees of freedom, rotation 
of the tool about its axis, provides the desired velocity of grinding (cutting) and is not 
related with the process for generation of £ a . 

The contents of this section cover the following topics: 

(1) Determination of equation of meshing between the tool surface E t and the generated 
surface E a . The equation of meshing provides the necessary condition of envelope 
existence to the family of surfaces. 
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(2) Determination of generated surface £ 5 as the envelope to the family of surfaces £ t 
swept out by the tool. Surface £ 5 coincides with the theoretical (ideal) surface £ p 
along the mean line L m and deviates from £ p out of L m . 

(3) Determination of deviations of £ ? from E P (in regions that differ from L m ) and mini- 
mizations of E s deviations for optimal approximation of E 5 to E p . 

(4) Determination of curvatures of E a that are required when the simulation of meshing 
and contact of two mating surfaces are considered. 

(5) Execution of required motions of E t with respect to E p by application of a multi-degree- 
freedom, computer numerically controlled machine. 

An effective approach for the derivation of the necessary condition of the envelope E 5 
existence is discussed. This method is based on the idea of motion of the Darboux-Frenet 
trihedron along L m , the chosen mean line of E p . 

An additional effective approach for determination of curvatures of generated surface E 3 
is discussed as well. This approach is based on the fact that the normal curvatures and 
surface torsions (geodesic torsions) of E fl are: (i) equal to the normal curvatures and surface 
torsions of E p along T m ; and (ii) equal to the normal curvatures and surface torsions of tool 
surface E t along the characteristic L g (the instantaneous line of tangency of E t and E s ). 

Mean Line on Ideal Surface E p 

The ideal surface E p is considered as a regular one and is represented as 

r» («,, K) 6 C 2 , («,,#„) 6 E (20) 

where (u p , 9 P ) are the Gaussian coordinates of £ p . 

The unit normal to £ p is represented as 
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N 0 


n p = 


No = 


dr v dr v 


|N P | ’ p du p ae p 


( 21 ) 


The determination of mean line on L m is based on the following procedure: 


(i) Initially, we determine numerically n points on surface E p that will belong approxi- 
mately to the desired mean line L m . 

(ii) Then, we can derive a polynomial function 


u pii 6 P>) = E j) ’ = 

i=i 


( 22 ) 


that will relate surface parameters (u p ,0 p ) for the n points of the mean line on E p . 

The mean line L m , tangent T p and unit tangent t p to the mean line are represented as 
follows 


^p) 1 



dr p du p 
du p dd p 



(23) 


The constraint for t p is that it must be of the same sign and differ from zero at the same 
intervals of interpolation. 

Tool Surface 

The tool surface E t is represented in coordinate system S t rigidly connected to the tool 
by the following equations 

x t = x t (u t )cos6 t , y t = x t (u t )sin0 t , z t = z t (u t ) (24) 
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The axial section of £ t obtained by taking 9 t = 0 represents a circular arc, or a straight 
line in the case when E t is a circular cone. Surface as shown in equations(24) of the tool is 
formed by rotation of the axial section of E t about the z t -axis. 

The surface unit normal is determined as 


n 4 


N « N = ^£i x ^i 

|N,r * dO t du t 


(25) 


Equation of Meshing Between £ t and £ 3 

Equation of meshing represents the necessary condition of existence of envelope £ 3 to 
the family of surfaces £ t that is swept out by the tool surface E t . 

The equation of meshing can be derived by using the equation 


Nf ) • v| ta) = 0 (26) 

Here: i indicates the coordinate system E; where the vectors of the scalar product are 
represented; is the normal to surface £ t ; is the relative velocity in the motion of 
E t with respect to £ 3 . 

Henceforth, we will consider two basic coordinate systems, S t and S p , that are rigidly 
connected to the tool surface E t and the ideal surface E p . In addition to £*, we will consider 
two trihedrons: 5(,(tf,, d fc , n b ) and Sy(ty, dy, ny). Trihedron is rigidly connected to £ f and 
coordinate system S t (figure 16). Here: Ob is the point of the chosen generatrix of £ f where 
the trihedron is located; t*, is the tangent to the generatrix at Ob’, n*, is the surface unit 
normal of £ t at Ob ; d& = x t&; vectors tj and d& form the tangent plane to E t at Ob- 
Trihedron 5/ moves along the mean line L m (figure 17); t f is the tangent to the mean line L m 
at current point M (figure 17); ny is the surface unit normal to E p at point M;dy = nyxty; 
vectors ty and dy form the tangent plane to £ p at point M. 

The tool with surface E t and trihedron S t moves along mean line L m of £ p and Ob 
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coincides with current point M of mean line L m . Surfaces £* and S p are in tangency at any 
current point M of mean line L m . The orientation of S with respect to St is determined 
with angle /3 that is varied in the process for generation. 

We start the derivations with the case when E t is a circular cone(figure 18). The angular 
velocity cj j of rotation of Sj with respect to S p is represented as 

= ( ttf — k n df + k g Tif)— (27) 


Here: t is the surface torsion (geodesic torsion), k n and k g are the normal and geodesic cur- 
vatures of surface E p at current point M of mean line L m , ds is the infinitesimal displacement 

along L m . 

The angular velocity j of trihedron S& is represented in S j as 


d0 \ , t d0 1 T ds 

n f = “f + s n i - l * k ‘ + * J it 


(28) 


The orientation of cone Et is determined by function 0(@p) and 

d0_ _ dP_ddp _ d0_ 1 / 2 g\ 

ds d9p ds d0 p |T p | 

where T p is the tangent to the mean line L m at current point M. 

The transformations of vector components in transition from S t to St, and St, to S/ are 
represented by 3 x 3 matrix operators Lf,t and L/t,. Here: 

cos 0 — sin0 0 

h/b = sin/3 cos 0 0 (^0) 

0 0 1 
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(31) 


Ut 


sin 7 1 cos 6 t sin 7 1 sin 6 t cos 7* 

sin 6 t — cos 0 

cos 7t cos cos 7 1 sin 9 t — sin 7* 


The cone surface E t is represented in S t as follows (fig. 18) 


r t = u t 


sin 7 1 cos 0 t sin 7 1 sin 0 t cos 7 1 


where (u t , 6 t ) are the surface parameters, 7 t is the cone apex angle. 
The unit normal to the cone surface is 


(32) 


n f 




cos 7 1 cos 0 t cos 7 t sin 6 t — sin 7 1 


The sought-for equation of meshing, necessary condition of existence of envelope 
represented in the form: 


(33) 
E s , is 


n 


(<) 

/ 



(34) 


where 


= L /t n t (35) 

The derivation of expression v ^ is simplified while taking into account the following 
considerations: 

(a) The relative velocity vector can be represented as 



= n 


(*)_(*) 
/ r / 


ds 

+ 7t t! 


(36) 
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Here, is the skew-symmetric matrix represented as 


Q 


(*) _ 


/ 


0 — u>3 LO 2 

CJ3 0 — 


— LO2 ^1 0 


Vector ft j is represented by 


ft j = f + d/ + — 


dtf 1 T ds 

t -K K+fe ] di 


( 37 ) 


(38) 


(b) Consider that point N on surface £ f (fig. 15) is the point of the characteristic (the 
line of tangency of £ £ and the generated surface £ 3 ). Certainly, the equation of meshing 
must be satisfied for point N. 

The position vector O/N can be represented as(figs. 15 and 18) 


0 f N = O t N - O t Of (39) 

Here, OjN is the position vector of point N that is drawn from the origin O t of S t to N; 
vector 0 t N is represented in St as 


O t N = u t e t = u t (sin 7 t cos 6 t i t + sin7 ( sin0 t jt + cos 7 t k £ ) 

where 



is the unit vector of cone generatrix O t N. 

Vector O^Oj (figure 18) is represented in S b as 


(40) 


( 41 ) 
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OtOf = It 4 


(42) 


where l t — O t O/|- 

Vector OjN is represented in 5/ using the matrix equation 


= u t L ft e t - ItLfbh 

(c) We represent now the equation of meshing as 


(43) 


n /' ’ v ;" ) = n ; !) ' [o w (u,L /t e, - W-jiM + (n f ■ t /) 


ds 

dt 


(44) 


(d) The further simplification of equation of meshing is based on the following rule for 
operations with skew-symmetric matrices [5]: 


A r B (a) A = C (,) 


(45) 


Here: B<*> and C<*> designate skew- symmetric matrices, A T is the transpose matrix for A. 
Considering that elements of B^ are represented in terms of components of vector 


b = [ b 2 b 3 j (46) 

we obtain that the elements of skew- symmetric matrix are represented in terms of 
components of vector c, where 


^T 


C\ C2 O3 


= A r 


i T 


bi hi 63 


(47) 


ds 

Using the above considerations and eliminating — , the final expression of equation of 

dt 

meshing can be represented as 
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(48) 


n 


W . = /(u t , 9 t , # P ) = titnf A^e t — 4- nfLj ( t/ — 0 


Here: 


A (,) — = 
dt 


L^ t f2 { / } L 


ft 


B<*>^ = L 
dt 


f2 ( , s) L fb 


fb**f 


(49) 


0 

-a 3 

a 2 


0 

-a i 

— a 2 

a l 

0 


(50) 


a i 
a 2 

a 3 


d p. 

t cos /3 sin 7 1 — sin /? sin 7t + (K + — ) cos j t 


t sin /? + k n cos 


df3. . 

t cos /? cos 7t — fc n sin /? cos 7 t — ( k g + -jj) sin j t 


(51) 



0 

—63 62 


61 


— £ cos fd + k n sin /? 

B (i) = 

^3 

0 -h 

5 

62 

r= — 

t sin (3 + k n cos /? 

_7 


— 6 2 

0 


&3 




(52) 


The family of characteristics L g and the instantaneous lines of tangency of £ ( and £ p are 
represented in St by the equations 


r f = r t (u t ,9 t ) , f(u t ,9 t ,9 p ) = 0 (53) 

where 9 p is the parameter of the family of L g . Taking 0 p = 9^ ( i = 1,2,..., n), we obtain 
the current characteristics on surface E ( . 
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It is easy to verify that the equation of meshing between E e and £ 3 is satisfied the current 
point M of mean line L m on the ideal surface £ p . This means that the characteristic L g 
intersects L m at point M, for which we can take 0 t = 0 since E t is a surface of revolution. In 
the case when £ t is a circular cone (figure 18), we can take for point M that u t = UtOb = U- 
The approach discussed above for the derivation of the equation of meshing can be easily 
extended for application in the more general case when the tool surface is a general surface 
of revolution. 

Determination of Generated Surface £ 5 

The ground surface E ff is generated as the envelope to the family of tool surface £ t ; 
surface £ s is represented in S p by the following equations 


r s P \ u p(Qp)iQpi u ti@t) — Lp/ r /^ + Tp \u p (6 p ),0p) , /(tX(, 0 t ,6 p ) — 0 (54) 

Here: f(u t ,0 t ,dp) = 0 is the equation of meshing; r^\u t ,9 t ) is the equation of tool surface 
E t represented in Sf, r^(up(& p ),O p ) is the vector function that represents in S p the mean 
line L m \ the 3x3 matrix operator L p / which transforms vectors in transition from Sj to S p 
is represented as 


tpx 

dp x 

n px 

k Vi = tpy 

dp y 

Upy 

r 

N 

dp Z 

Tlpz 


where 


(55) 
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n p = ± 


(57) 


dr p dr p 

de, 

dr p dr„ 

du p x de p 

d p = n p x t p ( 5 ®) 

The sign chosen in equation (57) must provide the direction of n p toward to the surface 
of the workpiece under consideration. 

Equations (54) represent in S p the generated surface E p in three-parametric form but 
with related parameters. Parameter u t is linear in equation of meshing when E t is a cone, 
therefore this parameter can be eliminated and the generated surface E 3 can be represented 

in S p as 


=r g = r g (9 p ,d t ) (59) 

Remember that surfaces S 3 and E p have a common line L m where they are m tangency. 
Surface E 3 is in tangency with E t along the instantaneous line L g that passes through current 
point M of line L m . The tangents to L g and L m lie in the plane that passes through M and 
is tangent to three surfaces (S P ,E 3 and E t ) simultaneously. 

Optimal Approximation of Generated Surface E 3 to Ideal Surface E p 

The procedure of optimal approximation of E s to E p is divided into the following stages: 
(i) design of grid on E p , the net of points, where the deviation of E s from S p will be 
determined; (ii) determination of initial function jd^\9 p ) for the first iteration; angle (3 
determines the orientation of the tool surface E t with respect to E p (figures 15 and 17); (m) 
determination of deviations of E s from E p with the initial function / 3 (1) (^ P ); 0 V ) optimal 
minimization of deviations. 

Grid on Surface E p . Figure 19(a) shows the grid on surface E p , the net of (n,m) points, 
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where the deviations of £ 3 from £ p are considered. The position vector is O p Qi t j = rj,* J ) 
(figure 19(b)). The computation is based on the following procedure: 

(i) The desired components Lij and Rij of the position vector r^’-d are considered as 
known. 

(ii) Taking into account that 

= fl? J = [xW)( Up .« p )] ! +[j,W)(u p ,«,,)] J (60) 

we will obtain the surface £ p parameters (u^'-d, 0('’-d) f or eac h grid point. 

Determination of Initial Function (3W (0 P ). The determination of j3^(0 p ) is based on 
the following idea: the instantaneous direction of tj (the tool generatrix) with respect to 
tangent tf to the mean line L m (figure 17) must provide the minimal value | 1 1 . Here: k^ 
is the relative normal curvature determined as 


ki r) = ki f) - ki p) 


(61) 


where k^ and k^ are the normal curvatures of surfaces E t and E p along tj,. In the case of 
nondevelopable ruled surface E p , vector tj, can be directed along the asymptote of E p . 

The requirement that | is minimal, enables to determine function Q^ l \6 p ) numerically. 

Since we need for further computations the derivative function ^ l \0 p ) is represented 

du p 

as a polynomial function that must satisfy the numerical data obtained for the chosen points 
of mean line L m . 

Determination of Deviations of E 3 from £ p . We are able at this stage of investigation to 
determine the equation of meshing between surfaces £ t and £ s , and surface S 5 as discussed 
above. The computation of deviations of E 5 from £ p at the grid points is based on the 
following considerations: 
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(i) Surfaces E p and E s are represented in the same coordinate system (S p ) by the following 
vector functions: 


r p (u p ,0p) , r g (9 g ,dt ) ( 62 ) 

(ii) The position vector and surface coordinates (itj* j) ,0£ j) ) are known for each point 
Q(*b) of the grid on surface E p . 

(iii) Point on surface E ? corresponds to point on surface E P . The surface E s 

parameters can be determined by using the following two equations 




(63) 


(iv) Due to deviations of E 3 from E p , we have that x£’ j) ^ 4* j) . The deviation of from 
E„ at the grid point is determined by the equation 


6i tj = nW-( r™-r™) (64) 

where nb ,J 4 is the unit normal to surface E p at the grid point Q p '^- 

The deviation can be positive or negative. We designate as positive such a deviation 


when Sij > 0 considering that nj/^ is directed into the “body” of surface E p . Positive devi- 
ations of E s with respect to E p provide that E s is inside of E p and surface E s is “crowned”. 


It is not excluded that initially the inequality Sij > 0 is not observed yet for all points 


of the grid. Positive deviations 8 t<J can be provided choosing the following options: 

(1) choosing a surface of revolution with a circular arc in the axial section instead of a circular 
cone; a proper radius of the circular arc must be determined. 

(2) changing parameter l t = |OtOfc| (figures 17 and 18); this means that the grinding cone 
will be displaced along tf, with respect to the mean line L m . 
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(3) varying the initially chosen function (3^(9 P ). 

Minimization of Deviations Consider that deviations ( i = j = 

of E ff with respect to E p have been determined at the (n, m) grid points. The minimization 
of deviations can be obtained by corrections of previously obtained function /?^( 0 P ). The 
correction of angle (3 is equivalent to the correction of the angle that is formed by the principal 
directions on surfaces £ t and E 5 . The correction of angle can be achieved by turning of 
the tool about the common normal to surfaces £ t and E p at their instantaneous point of 
tangency M*. 

The minimization of deviations j is based on the following procedure: 

Step 1: Consider the characteristic L g k, the line of contact between surfaces E t and £ s , 
that passes through current point Mk of mean line L m on surface E p (figure 20). Determine 
the deviations 8k between E t and E p along line L g k and find out the maximum deviations 
designated as and Points of L g k where the deviations are a maximum are 

designated as and N^. These points are determined in regions I and II of surface £ s 
with line L m as the border. The simultaneous consideration of the maximum deviations in 
both regions permits the minimization of the deviations for the whole surface E g . 

Note: The deviations of £ t from £ p along L g k are simultaneously the deviations of S 3 
from E p along L g k since L g k is the line of tangency of E t and £ s . 

Step 2: The minimization of deviations is accomplished by correction of angle (3k that is 
determined at point Mk (figure 20). The minimization of deviations is performed locally, 
for a piece k of surface £ fl with the characteristic L g k. The process of minimization is a 
computerized iterative process based on the following considerations: 

(i) The objective function is represented as 


F k 


= min ( 8< kL x 


+ 8 {2) 

' °kma 


with the constraint <5,j > 0. 


(65) 
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(ii) The variable of the object function is A0 k . Then, considering the angle 


1 = + Aft 


( 66 ) 


and using the equation of meshing with we can determine the new characteristic, the piece 
of envelope E^ and the new deviations. Iterations are required to provide the sought- for 
objective function. The final correction of angle f3 k we designate as (3 k ° pt \ 

Note 1: The new contact line (determined with /?£ 2) ) slightly differs from the real 

contact line since the derivative but not is used for determination Z$ . However, 

^gk very close to the real contact line. 

Step 3: The procedure discussed must be performed for the set of pieces of surfaces E 5 with 
the characteristic L gk for each surface piece. Remember that the deviations for the whole 
surface must satisfy the inequality S itj > 0. The procedure of optimization is illustrated with 
the flowchart shown in figure 21. 


Curvatures of Ground Surface E s 

The direct determination of curvatures of E s by using surface S 5 equations is a com- 
plicated problem. The solution to this problem can be substantially simplified using the 
following approach proposed by the authors: (i) the normal curvatures and surface torsions 
(geodesic torsions) of surfaces £ p and E s are equal along line L m , respectively; (ii) the nor- 
mal curvatures and surface torsions of surfaces E t and £ s are equal along line L g . This 
permits the derivation of four equations that represent the principal curvatures of surface 
2 in terms of normal curvatures and surface torsions of E p and E t . However, only three of 
these equations are independent (see below). 

Further derivations are based on the following equations: 


]c n = Jfc/ cos 2 q 4- ku sin 2 q = — (&/ 4- kn) 4- — (fcj — kn) cos 2 q 


(67) 
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t — 0.5(fcjj — kj) sin2<7 (68) 

Here: kj and kjj are the surface principal curvatures, angle q is formed by unit vectors 
e/ and e that is measured counterclockwise from e/ and e; e/ is the principal direction with 
principal curvature &/; e is the unit vector for the direction where the normal curvature is 
considered; t is the surface torsion for the direction represented by e. 

Equation (67) is known as the Euler equation. Equation (68) is known in the differential 
geometry as the Bonnet-German equation (see section 5). 

The determination of the principal curvatures and principal directions for E 3 is based on 
the following computational procedure (see section 5): 

Step 1: Determination of and for surface E 3 at the direction determined by the 
tangent to L m . 

The determination is based on equations (67) and (68) applied for surface E p . Recall 
that E p and E 3 have the same values of k^ and along the above mentioned direction. 
Step 2: Determination of k^ and . 

The designations k ^ and indicate the normal curvatures of E 3 and the surface torsion 
along the tangent to L g . Recall that k ^ and are the same for E t and E 3 along L g . We 
determine k^ and for surface E t using equations (67) and (68), respectively. 

Step 3: We consider at this stage of computation that for surface E 3 are known: and 

t^\ k ^ and t^ 2 \ for two directions with tangents r i and r 2 that form the known angle 
(fig. 22). Our goal is to determine angle qi (or q 2 ) for the principal direction e^ and the 
principal curvatures k\ 9 ^ and kff (figure 22). 

Using equations (67) and (68), we can prove that and ( i = 1,2) given for two 

directions represented by T\ and r 2 are related with the following equation 

*( 1 ) + *( 2 ) 

7m— i TiIT = cot '‘ < 69 > 

fi'Tl A'TX 
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Step 4: Using equations (67) and (68), we can derive the following three equations for 
determination of qi, k\ 9 ^ and k\j 


tan 2q\ ^ _ ^(i) cog 

( 70 ) 

k\ 9) = k^ - t (1) tan<fr 

(71) 

k\ 9 j = k^ + t (1) cot q x 

(72) 


Equation (70) provides two solutions for q x ( q J 2) = ^ + 90°) and both are correct. We 
choose the solution with the smaller value of q\. 

Numerical Example 1: Grinding of Archimedes’ Worm Surface 

The worm surface shown in figure 23 is a ruled undeveloped surface formed by the screw 
motion of straight line YW (\"KN\ = u p ). The screw motion is performed in coordinate 
system S p (figure 23(b)). The to be ground surface E p is represented in S p as 

r p = u p cos a cos 6 p i p + u p cos a sin 9 p j p + ( p6 p - u p sin a) k p (73) 

where u p and 9 P are the surface parameters. 

The surface unit normal is 


Thus: 



dr p dr p 
N p = X —— 
du p o9 p 


(74) 


33 



1 


n P /„.2 i ~.2 \0.5 


K + p 1 2 )°- ; 


p sin dp 4- Up sin a cos 0 P 
—p cos 6 p + Up sin or sin 9 p 


Up cos a 


As an example the following data will be used: 

Number of threads, N\, =2 


1 


(provided cosa^O) (75) 


Axial diametral pitch, P ax , = 8 (— ) 

m 

a = 20° 

The radius of the pitch cylinder = 1.125 (in) 

(i) The screw parameter is 


p = 


Ni 
2 P ai 


= 0.125 in. 


(ii) The lead angle is 


tan A„ = — = 


0.125 

1.25 


, A p = 5.7106° 


The mean line is determined as 


Fp(u m ,dp^ , u m 


r v + 


+ 


(•>-£) 


0.125 


r P - 


2 cos a 


COS Q 


= 1.3136 in. 


1 1.25 

where and — determine the addendum and dedendum of the worm. 

i ax * ar 

The worm is ground by a cone with the apex angle j t = 30°, and outside diameter 8 in.. 
The inside angle (3^ — —88.0121° provides the coincidence of both generatrices of the 
cone and the Archimedes’ worm. The maximal deviation of the ground surface from the 
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ideal surface E p with the above value of /3 (1) is 3 microns. 

The optimal angle /? (opt) = -94.6788° has been determined by the optimization method 
developed. The deviations of the ground surface from E p with the optimal /3 (opt) are 
positive and the maximal deviation has been reduced to 0.35 microns(figure 24). 

5. Condensed Information About Surface Curvature 

The contents of this section provide a condensed overview of about the basic equations 
of surface curvatures. For further explanation of the details please refer to the books by 
Nutbourne and Martin[6], Favard[7] and Litvin[2]. 

Osculating Plane 

Figure 25 shows spatial curve L X M L 2 . The osculating plane is the limiting position of 
such a plane that passes through curve points M\, M, and M 2 as M\ and M 2 approach M. 
The osculating plane for a curve at its regular point M is formed by the tangent to the 

curve and the acceleration vector for the same point. 

The osculating plane and the curve are in tangency of second order. The osculating plane 
is an exceptional tangent plane. The deviations of the curve from the osculating plane are 
of different signs on the two sides from the point of tangency, and the curve is above and 
below the plane (see points Li and L 2 in figure 25). An exception is the case when the point 
of tangency is a rectification point at which the second derivative r„ of a curve represented 
by r(s) is equal to zero. Here: s is the arc length of the curve. 

Space Curve and Surface Trihedron 

Henceforth, we will consider two trihedrons, the space curve trihedron and the surface 
trihedron. Each of the trihedrons is right-handed, formed by three mutually perpendicular 
vectors. The concept of space curve trihedron is discussed when a space curve is considered 
in the 3D space and the curve is not related to a surface. The concept of surface trihedron 
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and space curve trihedron are considered simultaneously when the space curve belongs to 
a certain surface and the curve inherits some of the properties of the surface to which it 
belongs. 

Space Curve Trihedron 

We consider a coordinate system that is rigidly connected to the curve. Position vector 
OC = r(s) determines the current point C of the curve (figure 25); s =MC is the length of 
the curve arc; M is the starting point. 

Consider that a small piece of curve L\MLi is located in the osculating plane IT 0 (figure 
25). Plane IUv is perpendicular to plane II 0 and passes through point M of the curve. 

We define the normal N to the curve as a vector that is perpendicular to the tangent 
to the curve. There is an infinite number of normals N to the curve at its point M . All of 
normals N belong to plane ITw since the unit tangent t is perpendicular to njv. For instance, 
vector N; is one of the set of curve normals (figure 25). Two normals of the set of normals 
must be specified: 

(i) the principal normal with the unit vector m that lies in the osculating plane IT 0 and is 
the line of intersection of planes IT 0 and II w (figure 25) and 

(ii) the binormal b that is perpendicular to t and m simultaneously. 

We may identify at each current point of the curve three mutually orthogonal vectors 
(figure 25): the tangent vector t, the principal normal m, and the binormal b. The ori- 
entation of these vectors in a fixed coordinate system is varied, depending on the location 
of the point on the curve. We may consider now a trihedron S c as a rigid body with three 
mutually perpendicular vectors e c (i c , j c , k c ) that form a right trihedron (figure 26). The 
origin of the trihedron moves along the curve, and the unit vectors i c , j c , k c represent t, m, 
b, respectively. Unit vectors t, m, b are taken at the current point of the curve where the 
origin of trihedron S c is located at this instant. 

The representation of unit vectors t, m, and b in terms of derivatives of vector function 
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r(s) is based on the following consideration: 

(i) Unit vectors t, m, and b form a right-hand trihedron (figures 25 and 26). Thus 

t = m x b , m = b x t, b = t x m (76) 

(ii) Unit vector t is directed along the tangent to the curve and therefore 


Vector r, is a unit vector since \ dr \= ds. 

(iii) The principal normal to the curve is perpendicular to the curve tangent t r 3 . The 

derivative r. 5 = — (r,) is perpendicular to r 3 , lies in the osculating plane and therefore the 

ds 

unit vector m of the principal normal is represented as 


m(s) = 


r„ 



(iv) Taking into account the expression for b in equations (76), we obtain the following 
equation for the binormal 


b(s) = t x m = 


r- x r. 


Frenet-Serret Equations 

The motion of the trihedren along a spatial curve can be represented in two components: 
(i) as a translational motion along the curve (the origin of the trihedron moves along the curve 
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and the unit vectors of the trihedron keep their original orientation), (ii) and as a rotational 
motion (the trihedron is rotated as a rigid body (to be coincided with the principal normal 
m c and the tangent t c to the curve at the curve neighboring point). 

Consider that the origin of curve trihedron coincides with point M of the curve and the 
unit vectors t c , m c and b c determine the instantaneous orientation of the trihedron(figure 
26). The neighboring point of the curve is N and \MN\ = ds , where s is the arc length of 
the curve. The unit vectors of the trihedron at N are determined as (t*, m* and b*), where 

t* = t c + t, e ds, m* = m c + m Sc ds, b* = b c + b 3C ds (77) 


Here: 


dt e dm c ^ db c 

= 1 ’ = ~3 ) bjc = ~z 

as as as 


(78) 


that are taken at point M. 

Frenet-Serret equations define t ic , m ic and b JC as follows (see References [6], [7] and [2]): 


t,c 
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o 
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m sc 

= 

rb c /C(jt c 
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o 
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-r m c 


0 -r 0 
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u 

s> 

1 


where k 0 and r are the curvature and torsion of the space curve at point M. It is evident 
that in the case of a planar curve, the unit vector b c is perpendicular to the plane where the 
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curve is located, b, c is equal to zero and the curve torsion t of a planar curve is equal to 
zero. 

Equations of k 0 and r for a Parametric Spatial Curve 

Consider that the spatial curve is represented by vector function r(6). After derivations 
we obtain (see References [2], [6], and [7]) 

r ge ■ m __ | r e x Tgg \ 

K ° " i r 2 , " |r e | 3 

(80) 

[(xgygg — xggye) 2 + (xgzgg — x gflZg) 2 + (t/gzgg — ygg^e) 2 ] 1 ^ 2 

(xff Tyf+ 2 |) 3/2 

The curvature k 0 obtained from equation (80) is always positive because the principal 
normal m c is located in the osculating plane and is directed to the center of curve curvature. 
The curve curvature k 0 can be also represented in the form 


k 0 = 



(81) 


Here: v r and a r are the velocity and acceleration of a point in its motion along the curve 
and are represented as follows 


dd 



(82) 


a r 


,<W * ,<P0. 

r " ( * ) +r,{ !? ) 


( 83 ) 


Obviously, the curvature k 0 can be also represented as 
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Tee • m 


(84) 


K 0 ~ 


The space curve torsion r is represented by the equation 


t(9) = 


(je x r gg ) • reee 
(re x Tee) 2 


(85) 


In the case of a planar curve, we have (r^ x Tee) • = 0 and r = 0. 

Surface Curve Trihedron 

Consider a regular surface £ that is represented by 

r(u, 9) e C 2 , r u xr # /0, (u,9) e A (86) 


A curve on £ is determined if in vector function r(u, 9) surface parameters are related 
with the equation 


f(u,9)= 0, f 2 + f 2 ^ 0 


(87) 


Figure 27 shows two curves, L n and Z 0 , that pass through the same surface point M and 
have the same tangent. Curve L n is a planar curve obtained by intersection of the surface 
by the surface normal plane that is drawn through the unit tangent t and the surface unit 
normal n. Curve L 0 is a spatial curve identified locally with the orientation of osculating 
plane, the curvature and the torsion of the curve. Considering that a spatial curve belongs 
to a surface, we may determine more parameters for the local identification of the curve. 


40 


We have introduced in above the curve trihedron ^(ic, jci ke) where i c t is the unit 
tangent, j c = m is the curve principal normal, and k c = b is the curve binormal (figures 
26 and 27(b)). In addition, we set up now the surface trihedron S/(i/, j y, k/) shown in 
figure 27(b). Here: iy = t is the unit tangent to the spatial curve, jy = d is the unit vector 
that is perpendicular to t and lies in the plane tangent to the surface at point M , ky — n is 
the surface unit normal. Subscript “/” indicates that the surface trihedron and its axes are 
considered. 

The unit tangent iy = i c = t is determined as 

t = -5- , T = r u + r e = r« - r 9 j- (/* ± 0) (88) 

| T | du Je 

The surface unit normal is represented as 

n = 7 ~~r , N = r u x r e (n = ky) (89) 


Changing the order in the cross product in equation (89), we can change the direction 
of n for the opposite one, and provide 8 < 90°, where 8 is formed by n and m. We remind 
that the direction of m is the same as r J4 (assuming that the curve is represented by r(s)) 
and cannot be chosen arbitrarily. Unit vectors t, d, and n form the right trihedron Sy, the 

surface trihedron. 

Bonnet-Kovalevski Equations 

Figure. 27(b) shows the curve and surface trihedrons whose common origin is located 
at the current point M of spatial curve L 0 . Consider now that the common origin of both 
trihedrons is moved along L a to the neighboring point N. Both trihedrons will keep the 
tangent t* to L 0 at point N as their common axis, but one of the trihedrons will be turned 
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with respect to the other one since the motion along L a will be accompanied with the change 
of angle 6 formed by vector m and n. Obviously, the unit vectors of the surface trihedron 
will change at N their orientation with respect to the orientation at M. Designating the unit 
vectors at N by t*, d’ and n*, we have 


t* = t(s) + t s ds, d“ = d(s) + d a ds, n* = n(s) + n s ds (90) 


where 


ts = d * = = ^( n ( 5 )) 


(91) 


Bonnet- Kovalevski equations express the derivatives t,, d 3 and n 3 in terms of tc g , rz n and 
t as follows (see References [6] and [2]). 


t s = K g d -f « n n = Kgjf + K n kf ] 


d, = -Kgt + tn = -Kgif + tkf > 
n s = — /c n t — td = —K n if — tjj 


(92) 


Here: K n , K g and t are the surface normal curvature, geodesic curvature, and the surface 
torsion, respectively. The concept of surface normal and principal curvatures is discussed in 
many books on differential geometry, but the determination and concept of K g and t requires 
additional explanation that is presented next in this report. 

Geodesic Curvature 

Frenet-Serret equations (92) yield that 
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t, = K 0 m 


(93) 


where *. is the curvature of a spatial curve; the curvature center lies in the osculating 
plane. Equations (92) and (93) yield that 


K„m = K,d + K„n 


(94) 


Equation (94) can be interpreted follows: 

(1) Figure 29 shows a spatial curve L , , on surface S. Unit vectors t, d, and n represent the 
surface trihedron (figures 29 and 27(b)). Here: t is the tangent to curve L 0 ; d lies in the 
tangent plane and is perpendicular to t; n is the surface unit normal. Unit vector m is the 

principal normal to L 0 and lies in the osculating plane. Vector r ja = K 0 m. 

(2) Consider now that the spatial curve L. is projected on the tangent plane T and normal 
plane N, respectively. The projections are designated by L T and L N . We emphasize that 
there is no difference between L„ (figure 29(b)) and L„ (figure 28) if they are considered 
locally. Both curves have the same normal curvature at the point of tangency M. 

(3) Vector ,c„m is represented as the sum of two vectors: a s d and *„n. The scalar k, 
represents the curvature of curve L T , and the scalar represents the curvature of curve I„. 

(4) Equation (94) yields two relations 


K 0 (m • n) = k 0 cos 8 — K, n 


(95) 


Kg = r„ ■ d = K 0 sin 8 


(96) 
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where S is the angle formed by vectors m and n that determines the orientation of osculating 
plane with respect to the normal plane. Equations (95) and (96) relate curvatures k 0 and 
K n and angle 6. 

The direct determination of geodesic curvature of a spatial curve represented on a surface 
is based on the following equations: 

(i) Consider that the surface is represented by the vector function r(u,0). 

(ii) A spatial curve is represented on the surface as 

r(u(0),0) (97) 


where u(6) is the known function. 

(iii) The tangent to the curve is represented as 


_ du 

T = r 'le + r " 


(iv) The unit normal to the surface is represented as 


N 

n = N = r„xr ( 


(v) An auxiliary parameter a is represented as 


, du . , _ du 

a = r uu (— ) +2r ue — + 


(98) 


(99) 


(100) 


(vi) The final expression for K g is 
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( 101 ) 


K 


9 


T-(a x n)- | N | 


<Pu 

J(P 


I T | 3 


Surface Torsion 

The surface torsion t can be represented by the equation 

( = T + 6, = T + ~(S) (102) 

as 

Thus, the surface torsion depends on the torsion r of the spatial curve along which the 
origin of two trihedrons is moved, and on the derivative S s where 6 is the angle formed by 
the unit vectors m and n of the trihedrons. 

The geometric interpretation of the surface torsion may be based on the concept of the 
geodesic line (see References [6] and [2]). A spatial line on the surface is the geodesic one if 
the principal normal m at any curve point M coincides with the surface normal at M. The 
geodesic curvature of the geodesic line at any curve point is equal to zero. 

It was proven in differential geometry that the surface torsion t is the curve torsion of 

the geodesic line. 

A simple method for computation of the surface torsion is based on the equation that 
has been proposed by Sophia Germain and Bonnet (see References [6], [2] and [(])• This 
equation is 


t = 0.5(/c/j - /cj)sin2 q 


(103) 


Here: kj and k u are the principal curvatures of the surface at point M on the principal 
directions with the unit vectors ej and e//(fig. 30); q is the angle formed by e/ and t. 
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Using equation (103) and Euler’s equations that relate the principal curvatures and nor- 
mal curvatures, we may determine the solutions to the following two problems (Ref. [2]): 
Problem 1: Consider that two directions in the tangent plane determined with unit vectors 
t*0 and tU) are given (figure 31). The angle fi formed by U 1 ) and U 2 ) is known. Also the 
following are known: (i) the normal curvatures /c^ and for directions of U 1 ) and U 2 \ 
and (ii) the surface torsion t' 1 ) given for direction U 1 ). 

The goal is to determine the principal curvatures kj and kjj for directions of U 1 ) and 
t* 2 \ and angle q ^ (or ? (2) ). 

The solution to this problem is represented by the following equations [2] 


tan2</ 1) = 

t^(l — cos2/z) 

(104) 

Kn ^ — /cl 1 * — tU) sin 2fi 

Kl = 

/c^-t^tangW 

(105) 

Kn = 

/cW-M^cotgW 

(106) 


Problem 2: Consider as given t^, U 2 ) (figure 31), and k^. The goal is to relate the 
surface torsions for directions of t^ and t^. 

The sought-for relation is represented by the equation [2] 


*(i) +*(2) 

15) (ij--cot/x = 0 (107) 

Kn Kn 


where fi is the angle formed by t^ and U 2 ). 

Numerical Example 2: Determination of the geodesic curvature K g and surface 
torsion t of Archimede’s worm surface 
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Archimede’s worm surface is represented as 

u cos a cos 9 

u cos or sin 6 (108) 

p9 — u sin a 

Here u and 6 are the Gaussian coordinates, a is the pressure angle and p is the screw 
parameter. A helix on the Archimede’s worm surface is a spatial curve obtained by inter- 
section of the worm surface by a cylinder of radius r t . Our goal is to determine the geodesic 
curvature K g of the helix and the surface torsion of the Archemede’s worm surface. 

(1) Geodesic curvature K g 

Taking into account the r\ + r 2 y = rf, we can represent the helix as follows 

[ r ,• cos 6 1 

u = r(9) 

cos a 

| p9 — r{ tan a J 

The tangent to the helix can be represented as 

T = — = (— r,sin0 r,-cos0 p) T (HO) 

o9 

The unit surface normal is represented by equations (99) and (75) 
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n = 


y/p + U 2 


p sin 9 + u sin a cos 0 
— p cos 0 -f u sin a sin 9 


u cos or 


( 111 ) 


The auxiliary vector a is (see equation (100)) 

a = (— ri cos 9 — r^sin# 0) r 


( 112 ) 


Equations (101) and (110) to (112) yield the following expression for the geodesic curva- 


ture 


T • (a x n) 


r,- cos a 


\/(p*~+"r?)(p*c os^aTrfj 


(113) 


(2) Surface torsion t 

From Reference [2] (F. L. Litvin, 1993), the principal curvatures and principal directions 
at a surface point can be represented by the following equations: 


Ki = 


e t = 


Lhi + M 
Ehi + F ’ 

(• = i, a) 

(114) 

r u hi + r e 
|r u /i,- + r e | 

(* = i, in 

(115) 


The coefficients and the partial derivative in the case of the Archimedes’ worm can be 
expressed as follows: 
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1 = 0 , M = -TT^T\ X’ F=-psina, E = 1 

(u 2 + p 2 ) 2 

_ -u 2 sin q ± (u 4 sin 2 a + 4p 2 u 2 + 4p 4 )^ .. _ j 

li = 2p ’ * ’ ' 


(116) 


(117) 


r u = (cos a cos 6 cos a sin 5 — sin a) 1 


(118) 


1 


re = 


y/p* + u 2 cos 2 a 


:( — xi cos a sin 0 u cos a cos 6 p) 1 


(119) 


Angle q that is formed by tangent T and ej is 


_ x ,ej-T 
q = cos (-j^p) 


(120) 


Considering that «/, K It and angle q are given, we can obtain the surface torsion along 
the tangent T as 


t = 0.5(/cjj — K/)sin2 q 


( 121 ) 


(3) Computation results 

The to be computed point is located on the helix that belongs to the pitch cylinder of 
the worm. The z-coordinate of the helix point is equal to zero, and the Gaussian coordinates 

are 


= — , 0 = — tana 

cos a p 
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The design parameters are the same as in Numerical Example 1, i.e. 


r p = 1.25m, p = 0.125m, a = 20° 

The results of computation are 

k„ = 0.8257—, t = —0.1540—, ki = —0.3283—, «// = 0.0227— 

3 in in m m 

6. Conclusion 

From the analytical study presented in this report the following conclusions can be drawn: 

(1) The kinematics of two CNC machines with 6 degrees- of-freedom has been described. 

(2) The preliminary results of investigation by TCA of the sensitivity of helical gears and 
worm-gear drives to misalignment are represented. 

(3) A new method for grinding of a gear tooth surface with optimal approximation to the 
given surface is proposed. 

(4) An algorithm for the execution of motions of a CNC machine for the surface generation 
has been developed. 
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Fig. 4 Transmission error of helical gears: error of pinion 
lead angle on pitch cylinder is AA p i = 5.0arc - min 
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Fig. 6 Shift of bearing contact due to change of center 
distance A E and crossing angle At 
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Fig. 9 Schematic of “Phoenix” machine 
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Fig. 10 Coordinate systems applied to “Phoenix” machine 
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Fig. 11 “Star” CNC machine 
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i - tilt angle 


Fig. 12 Pinion head-cutter 
63 







Fig. 15 Installment and orientation of tool surface £* 
with respect to ideal surface £ p 
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Fig. 16 Tool surface 
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Fig. 17 Orientation of trihedron Sb with respect to 5/ 
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Fig. 19 Grid on surface E p 




Fig. 20 Determination of maximum derivations 
along line L g * 
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Fig. 21 Flow chart for optimization 
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Fig. 22 To determination of principal directions 
of generated surface 
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Fig. 23 Surface of Archimedes’ worm 
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Fig. 24 Deviations of the ground surface from 
ideal surface S p of Archimedes’ worm 
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Fig. 25 Osculating and normal planes 
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Fig. 26 Curve trihedron 
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Fig. 27 Surface trihedron 
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Fig. 28 Surface normal section and surface spatial curves 


79 


Tangent plane T 









REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions searching existing data sources 
gathering and maintaining the data needed, and completing and reviewing the collection ol information. Send comments regarding this burden estimate or any other aspect of this 
osllection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704*0188), Washington, DC 20503. 


1. AGENCY USE ONLY {Leavebfank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

April 1994 Final Contractor Report 


4. TITLE AND SUBTITLE 


Generation of Gear Tooth Surfaces by Application of CNC Machines 


6. AUTHORS) 


F.L. Litvin and N.X. Chen 


5. FUNDING NUMBERS 


WU-505-62-36 

1L162211A47A 


7. PERFORMING ORGANIZATION NAM E(S) AND A DDR ESS(ES) 

University of Illinois at Chicago 
Chicago, Illinois 60680 


a PERFORMING ORGANIZATION 
REPORT NUMBER 


E-8712 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

Vehicle Propulsion Directorate 
U.S. Army Research Laboratory 
Cleveland, Ohio 44135-3191 
and 

NASA Lewis Research Center 
Cleveland, Ohio 44135-3191 


11. SUPPLEMENTARY NOTES 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR- 195306 
ARL-CR-145 


Project Manager, Robert F. Handschuh, Propulsion System Division, organization code 2730, (216) 433-3969. 


12a. DISTRIBUTION/AVAILABIUTY STATEMENT 

Unclassified - Unlimited 
Subject Category 37 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

This study will demonstrate the importance of application of CNC machines in generation of gear tooth surfaces with new 
topology. This topology decreases gear vibration and will extend the gear capacity and service life. A preliminary investi- 
gation by a tooth contact analysis(TCA) program has shown that gear tooth surfaces in line contact (for instance, involute 
helical gears with parallel axes, worm-gear drives with cylindrical worms etc.) are very sensitive to angular errors of 
misalignment that cause edge contact and an unfavorable shape of transmission errors and vibration. The new topology of 
gear tooth surfaces is based on the localization of bearing contact, and the synthesis of a predesigned parabolic function of 
transmission errors that is able to absorb a piecewise linear function of transmission errors caused by gear misalignment. 
The report will describe the following topics: (1) Description of kinematics of CNC machines with 6 degrees-of-freedom 
that can be applied for generation of gear tooth surfaces with new topology. (2) A new method for grinding of gear tooth 
surfaces by a cone surface or surface of revolution based on application of CNC machines. This method provides an 
optimal approximation of the ground surface to the given one. This method is especially beneficial when undeveloped 
ruled surfaces are to be ground. (3) Execution of motions of the CNC machine. The solution to this problem can be 
applied as well for the transfer of machine-tool settings from a conventional generator to the CNC machine. The devel- 
oped theory required the derivation of a modified equation of meshing based on application of the concept of space 
curves, space curves represented on surfaces, geodesic curvature, surface torsion, etc. Condensed information on these 
topics of differential geometry is provided as well. 


14. SUBJECT TERMS 

Gears; Gear teeth 

IS. NUMBER OF PAGES 

84 

16. PRICE CODE 

A05 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

20. LIMITATION OF ABSTRACT 

OF REPORT 

OFTHISPAGE 

OF ABSTRACT 


Unclassified 

Unclassified 




NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 
















National Aeronautics and 
Space Administration 

Lewis Research Center 
21000 Brookpark Rd. 
Cleveland, OH 44135-3191 



